13 research outputs found

    Accurate Profiling of Microbial Communities from Massively Parallel Sequencing using Convex Optimization

    Full text link
    We describe the Microbial Community Reconstruction ({\bf MCR}) Problem, which is fundamental for microbiome analysis. In this problem, the goal is to reconstruct the identity and frequency of species comprising a microbial community, using short sequence reads from Massively Parallel Sequencing (MPS) data obtained for specified genomic regions. We formulate the problem mathematically as a convex optimization problem and provide sufficient conditions for identifiability, namely the ability to reconstruct species identity and frequency correctly when the data size (number of reads) grows to infinity. We discuss different metrics for assessing the quality of the reconstructed solution, including a novel phylogenetically-aware metric based on the Mahalanobis distance, and give upper-bounds on the reconstruction error for a finite number of reads under different metrics. We propose a scalable divide-and-conquer algorithm for the problem using convex optimization, which enables us to handle large problems (with 106\sim10^6 species). We show using numerical simulations that for realistic scenarios, where the microbial communities are sparse, our algorithm gives solutions with high accuracy, both in terms of obtaining accurate frequency, and in terms of species phylogenetic resolution.Comment: To appear in SPIRE 1

    Use of secondary forests by understory birds in a fragmented landscape in central Amazonia

    Get PDF
    Rates of deforestation in the Brazilian Amazon have increased since 1991 and forecasts are not optimistic about the slowing of this process. Some authors believe that the Amazon may be experiencing a massive process of species extinction. However, the deforestation is accompanied by the expansion of secondary forests that are established in the abandoned areas. The trend is an increase in secondary forests cover, resulting in a mosaic of primary forest (FP) and fragments separated by an array of secondary forests (FS). In this scenario, the prediction of a massive extinction could be wrong if many species could survive in the secondary forests. To assess the importance of FS for the understory birds we sampled areas in regeneration and a continuous forest of a fragmented landscape. We conducted mist netting (24 nets/day) for six consecutive days/month, for 8 months (May-November) in 2009. Some forest species as do not seem to be adapted to the secondary forest environment and their occurrences are restricted to continuous forest environments. But most focal species showed no significant difference in apparent survival rates between the enviroments, suggesting that these species inhabit the secondary forest and the primary forest similarly. Because most of the matrix in fragmented landscapes are composed by secondary forests, such results highlights the conservation value that these habitats present in the long term. Thus, FS should be regarded as dynamic matrix that not only allows the movement of individuals but also function as habitat for many species typical of FP.Na Amazônia, as taxas de desmatamento crescem desde 1991 e as previsões não são otimistas quanto à desaceleração desse processo. A devastação da floresta é acompanhada de uma expansão de florestas secundárias (FS) que se estabelecem nas áreas abandonadas. A tendência é um aumento de florestas secundárias, resultando num mosaico de floresta contínua e fragmentos separados por uma matriz de FS. Nesse cenário, autores acreditam que a Amazônia pode passar por um processo massivo de extinção de espécies. Por outro lado, a previsão de um processo massivo de extinção pode ser equivocada, pois muitas espécies florestais poderiam sobreviver nas florestas secundárias. Para avaliar o valor das florestas secundárias para espécies florestais amostramos por oito meses com redes de neblina uma capoeira (FS) em regeneração e uma floresta primária (FP) de uma paisagem fragmentada. Algumas espécies não foram capturadas na capoeira e aparentemente evitam esse tipo de hábitat. No entanto, a maioria das espécies do grupo focal não apresentou diferença na sobrevivência aparente entre os ambientes, o que nos indica que estão habitando a capoeira e a floresta primária da mesma forma. Na realidade amazônica, onde grande parte da matriz é composta por floresta secundária, a matriz tem valor para conservação e deve ser analisada como um elemento dinâmico que não apenas permite a movimentação de indivíduos, mas também serve de hábitat para muitas espécies de floresta primária. Mas ressaltamos que é fundamental a preservação de áreas de floresta primária que servirão de fonte às florestas secundárias adjacentes
    corecore